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Abstract 

The quark-gluon to hadron transition affects the evolution of cosmological 
perturbations. If the phase transition is first order, the sound speed vanishes 
during the transition, and density perturbations fall freely. This distorts 
the primordial Harrison-Zel'dovich spectrum of density fluctuations below the 
Hubble scale at the transition. Peaks are produced, which grow at most 
linearly in wavenumber, both for the hadron-photon-lepton fluid and for cold 
dark matter. For cold dark matter which is kinetically decoupled well before 
the QCD transition clumps of masses below 10~ 10 Mq are produced. 
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QCD makes a transition from a quark-gluon plasma at high temperatures to a hadron gas 
at low temperatures. Lattice QCD simulations give a transition temperature T* ~ 150 MeV 
and indicate a first-order phase transition for the physical values of the u,d,s-quark masses [|TJ. 
The relevance of the QCD transition for cosmology, especially for big-bang nucleosynthesis 
0, has been discussed before, but the focus was on effects of bubble formation 0|§. In 
this paper we look at matter averaged over scales A much larger than the bubble separation. 
We show that for a first-order phase transition the sound speed c s = (dp/dpYJ 2 drops to 
zero for these wavelengths when the transition temperature T* is reached, stays zero for the 
entire time until the phase transition is completed, and afterwards suddenly rises back to 
c s « c/ v^- In contrast the pressure stays positive and varies continuously, although it goes 
below the radiation fluid value p = p/3. Since c s is zero during the transition, there are no 
pressure perturbations, no pressure gradients, no restoring forces. Pre-existing cosmological 
perturbations, generated by inflation || with a Harrison- Zel'dovich spectrum ||, go into 
free fall. The superhorizon modes (at the time of the transition) remain unaffected. The 
subhorizon modes develop peaks in Sp/p which grow with wavenumber k > k+, where k^ hys ~ 
Hubble rate H at the end of the QCD transition. The details of this growth depend on the 
QCD equation of state near T+. We analyze two cases: First we use the bag model 0, 
which gives a simple parameterization and allows a simple discussion of the effects. It gives 
a maximal latent heat, and produces peaks in 5p/p which grow linearly in k. Next we 
use lattice QCD results fl||^§], which indicate a smaller latent heat, and we fit s/T 3 = 
C\ + C^l — T^/T) 1 / 3 above T*. This produces peaks in Sp/p which grow as k 3 / A . 

The sound speed (for wavelength A much larger than the bubble separation), c s = 
(dp/dp) l J 2 , must be zero during a first-order phase transition of a fluid with negligible 
chemical potential, since the fluid must obey 

» + P=T%, (1) 

according to the second law of thermodynamics. Because the energy density p is discontinu- 
ous in temperature at T* for a first-order phase transition, the pressure p must be continuous 
with a discontinuous slope. As the universe expands at fixed temperature T* during the phase 
transition, p as a function of time slowly decreases from p+(T*) to p_(T*), p stays constant 
at p{T+), and therefore c s is zero during the whole time of the phase transition. 

The interaction rates in the QCD-photon-lepton fluid are much larger than the Hubble 
rate, T/H 1, therefore we are very close to thermal and chemical equilibrium, the QCD 
transition is very close to a reversible thermodynamic transformation, and the entropy in 
a comoving volume is approximately conserved. Estimates show that supercooling, hence 
entropy production, is negligible, (T* — T supercoo i ing )/T + ~ 10~ 3 [|10j. Bubble formation during 
the QCD phase transition is unimportant for our analysis, estimates give a bubble separation 
£b ~ 1 cm while the Hubble radius at the QCD transition is Rh ~ 10 km, therefore 
ibjRu ~ 10~ 6 . We shall analyze perturbations with A 3> £&. 

In the bag model J7J it is assumed that for T > the quark-gluon plasma (QGP) obeys 

P QGP (T) = P^l(T) - B , (2) 

where Pqqp{T) = (7t 2 /90)(7q GP T 4 , g* is the effective number of relativistic helicity states, 
and B is the bag constant. We include u,d-quarks and gluons in the quark-gluon plasma, 
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7, e, p, and 3 neutrinos in the photon-lepton fluid (7L), and for T < we have a hadron gas 
(HG) of pions. We treat the pions as massless and ideal, because their contribution is small 
anyway, S'hg/S'qgp = 3/37 and S'hg/^l = 3/14.25. p follows from Eq. (Q) via the second law 
of thermodynamics, Eq. ([I]), and s from s = dp/dT. This gives for the quark-gluon plasma 

PQGp(T)=pg£(T)+B (3) 
*QGp(T)=agSSp(r) . (4) 

The bag constant is determined by the critical temperature T* via pqqp(T + ) = Phg(^*)- 

The latent heat, L = T+As, should be compared with the difference in entropy between 
an ideal HG and an ideal QGP. This defines the ratio Rl = L/(T*As) ldeal . The bag model 
gives Rl = 1. Lattice QCD indicates that the transition is first order both for quenched 
QCD (no dynamical quarks) and for QCD with three quarks and physical masses |L]. For 
the latter case neither the value of the latent heat nor the equation of state are available. 
Quenched QCD gives L/T* « 1.4 0, which implies R L « 0.2. We fit the shape of the QCD 
entropy to quenched QCD data [§] for T > by 



QGP = i + 

„ideal 
6 QGP 



T x7 
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T ' J fi'QGP 



A ' ;< [I-Rl) , (5) 



where Ag* = gq GP — g^ G and Rl = 0.2. A good fit for our purpose is obtained for 7 6 
(0.3,0.4). We fix 7 = 1/3. 

The growth of the scale factor during the c 2 s = part of the QCD transition, a + /a_, 
follows from the conservation of entropy in a comoving volume, 



a + 



Ag" 

i + R L - 



Salter 




Rl — 1 

i? L = 0.2 ' 1 j 



taking into account photons, leptons, and hadrons in g* fter . Fig. |l| shows the evolution of the 
sound speed with the scale factor a. Above T*, the sound speed in the bag model has the 
value for an ultrarelativistic ideal gas, c s — 1/ v^3, because the bag constant drops out when 
forming dp/dT and dp/dT in Eqs. (H) and @. The sound speed vanishes for about a third 
of a Hubble time for Rl = 1 and for a tenth of a Hubble time if Rl = 0.2. The pressure 
does not drop all the way to zero, it drops to p^iT*) +Phg(^*)- 

The evolution of linear cosmological perturbations through the QCD transition is ana- 
lyzed in the longitudinal sector (density perturbations) for perfect fluids. We choose a slicing 
E of space-time with unperturbed mean extrinsic curvature, 5[tri£y(E)] = 0. This implies 
that our fundamental observers, which are defined to be at rest on the slice E, w(obs) = n(E), 
have relative velocities, which in the mean over all directions follow an unperturbed Hubble 
flow. If the coordinate choice (gauge choice) is such that the time coordinate t is constant 



on the slices E, the gauge is fixed to be the uniform expansion (Hubble) gauge |11]. As fun- 
damental evolution equations for each fluid we have V^T^ = 0, i.e. the continuity equation 
and (in the longitudinal sector) the 3-divergence of the Euler equation of general relativity, 

d t e = -3H(e + 7r)- A<0-3#(p + p)a (7) 
d t ip = -3#V - n-{p + p)a , (8) 
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where e = 5p, it = dp, p = p , p = p , Vip = S — momentum density (Poynting vector), a = 
lapse function. The system of dynamical equations is closed by Einstein's -Rgg-equation, the 
general relativistic version of Poisson's equation, 

(A + 3H)a = 47rG(e + 3tt) , (9) 

together with the equation of state. Equations (|7|) - (0) define our general relativistic Cauchy 
problem in linear perturbation theory in the longitudinal sector with initial data (e, ip) freely 
chosen on Ej. These three equations are the Jeans equations extended to general relativity 
in the longitudinal sector. In all three of them the mean over all directions is taken. This 
fact matches our slicing condition that in the mean over all directions the relative velocity 
of our fundamental observers is unperturbed. Therefore the uniform expansion (Hubble) 
gauge could be called the 'longitudinal Jeans gauge'. 

It is convenient to work with the dimensionless variables 5 = e/p (density contrast), 
if; = k phjs ip/p (~ peculiar velocity) and with conformal time, ()' = d v = adt. In our 
numerical analysis we have used the exact general relativistic equations, but it is instructive 
to look at the subhorizon approximation, A phys <C H~ x , where one can drop H in the general 
relativistic Poisson equation ([5]) and the time dilation term (last term) in the continuity 
equation (^|). Furthermore, if we take the limit in which the QCD transition time is much 
shorter than the Hubble time, (£+ — £_) <C H , and if we integrate the equations during 
a correspondingly short time interval, we can drop the remaining terms proportional to H, 
and Eqs. (0) - (||) simplify to 

5' = kip 

= - c 2 s k8- {l+w)ka (10) 

(^) 2 *=-^ +3 ^' 

where w = p/p. With these approximations the evolution of a mode k of cosmological 
perturbations can be solved analytically. 

The origin of large peaks in 8p/p for k ^> k*, where k phjs ~ H at the transition, is 
easily understood in the bag model. For the dynamics of the radiation fluid (QCD, photons, 
leptons) one can neglect cold dark matter (CDM) since Pcdm/Prad = a / a equaiit y ~ 1CT 8 . 
The radiation fluid in each mode makes standing acoustic oscillations before and after the 
QCD transition with gravity negligible and with equal amplitudes of 5 and s/Zip. The 
solution before the QCD transition is 8(r]) = Ai n cos[u (rj — r^_) — and ip = S'/k, where 
to = kc s = kj a/3, and 95 _ is the phase of the acoustic oscillation at r^_, i.e. at the beginning 
of the QCD transition. During the QCD transition the sound speed is zero, there are no 
restoring forces from pressure gradients, the radiation fluid goes into free fall. But during 
this free fall gravity is again negligible for the radiation fluid, if (t+ — £_) H -1 . This is 
inertial motion in the sense of Newton. The peculiar velocity is constant in time, ip{v) = 
and the density contrast grows linearly in time with a slope k, S(rj) = <5_ + k(r] — 7]J)ip_. 
Thus, the final amplitude A + grows linearly in k modulated by sin(y9_), which produces 
peaks in the spectrum. The height of these peaks is 

(V+ ~ V-)k _ k 

/q ~ h 1 ^ ' 

peaks V "J "'l 
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for k ^> k\. The usual free fall growing mode at subhorizon scales behaves totally different: 
ip/8 ~ H/k phjs <C 1 and 5 + /5- = (i] + /i]_) 2 m 2. In our case the initial peculiar velocity ip 
from the acoustic oscillations of the radiation fluid is enormously larger than in the usual free 
fall growing mode for k 3> k*. COBE observations |12[ normalize the subhorizon spectrum 



of density perturbations as {Sp/p)\ ~ 10 -4 , if there is no tilt in the spectrum. Hence in 
the bag model radiation-fluid modes with k/k\ ^ 10 4 go nonlinear by the end of the QCD 
transition. 

Our numerical results for the spectrum of density perturbations from the lattice QCD 
fit Eq. (W) are given in Fig. [2[ We show the enhancement of the amplitude A R ad = (^rad + 
S^rad) °f the acoustic oscillations of the radiation fluid after the transition compared 
to the amplitude without transition. For CDM we show the amplitude Acdm = |^cdm| 
at T*/10 compared to Arad without transition. In both cases we obtain peaks over the 
Harrison-Zel'dovich spectrum of primordial adiabatic density fluctuations. The modes k 
(horizontal axis) are labeled by the CDM mass contained in a sphere of radius A/2 = %/k. 
The positions of the first few peaks and dips are the same for the lattice QCD fit and for 
the bag model. In both cases the beginning of the peak-dip structure is at fci, see Eq. flTTD , 
which is ~ k^. k\ corresponds to Mp DM « 9 x 1O~ 9 M . The peaks grow as (k/k 2 ) 3 ^ for 
k ^> k 2 , see Eq. (O). The radiation energy inside Ai/2 is ~ 1M Q , but it gets redshifted as 
^RAD( a ) ~ («equaiit y /a)^CDM- The time evolution of a subhorizon mode is shown in Fig. |3|. 
During c s = 5rad grows linearly as in the bag model. Above T* the evolution differs 
because c s < l/v3. 

The slower rise in k for the lattice QCD fit can be understood by a WKB analysis. With 
the same approximations as above, Eqs. ( |T0| ) reduce to 

5" + c 2 s k 2 5 = . (12) 

Under the WKB condition, \dc s /dr]\/c s <C u = c s k, the solution reads S = 
An(3c^)~ 1//4 cos(/c / c s drj). Note that the peculiar velocity (~ if)) decreases with c l J 2 , which 
can be seen in Fig. |3|. Just above T* the sound speed may be approximated by c s oc (rj- —rf). 
This gives the solutions of Eq. ([12]) as 5 oc z 1 / 4 J±iu(z) with z = c s k(r)- —rj)/2. The normal- 
ization for large k is provided by the WKB solution. At r]_ one matches the linearly growing 
solution in the regime of vanishing sound speed (r/_ < i] < The final amplification is 



A, 



71 



peaks Z,J - \ 4/ 



3k 




( V+ - V -)k=[-\ , (13) 



for k ^> k 2 . k 2 depends on the duration of the c s = regime and thus depends on R L . 
For Rl = 0.2 it corresponds to a mass M 2 = 2 x 1O _1O M , which is indicated in Fig. |2| 
together with the asymptotic envelope (k/k 2 ) 3 ^. Without tilt COBE normalized modes 
with k/k 2 > 10 5 go nonlinear by the end of the QCD transition. The perturbations in the 
radiation fluid will get wiped out by collisional damping from neutrinos at temperatures 
below T+ but above 1 MeV. 

For cold dark matter we consider any non-relativistic matter which decouples kinetically 
well before the QCD transition. The neutralino, most likely the lightest supersymmetric 
particle [O, is weakly interacting. Thus it decouples kinetically around T ~ 1 MeV and is 
excluded to make our CDM. It would belong to the radiation fluid at T*. Candidates for 
our CDM are axions or primordial black holes. 
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CDM falls into the gravity wells generated during the transition by the radiation fluid. In 
the bag model this leads to peaks in CDM which grow linearly for k 3> k%, i.e. ^cdm^cdm _ 
[H(t+ - t_)/2] 2 5* AD . After the transition S CDM grows logarithmically. For our lattice QCD 
fit the CDM peaks are shown in Fig. |2j. An analytic analysis of the CDM evolution will be 
presented in a longer paper. 

The implications of these peaks above the Harrison-Zel'dovich spectrum generated in a 
standard scenario with a first order QCD transition are: 

1) For CDM which is kinetically decoupled well before the QCD transition clumps with 
^cdm ~ 1O _1O M are produced. They go nonlinear after equality and virialize by violent 
gravitational relaxation. Assuming a COBE normalized spectrum with tilt n — 1 = 0(0.2) 
and 3a peaks the size of 10~ 10 M o clumps is ps 14 AU (1 AU). 

2) Big-Bang Nucleosynthesis will not be affected by nonlinear acoustic oscillations of the 
radiation fluid for k/ki > 10 4 (bag model) resp. 10 6 (lattice QCD fit), because they are 
wiped out by collisional damping from neutrinos before Big-Bang Nucleosynthesis. 

3) Primordial black hole formation is unlikely, because the nonlinear acoustic oscillations in 
the radiation fluid are generated far below the Hubble scale at the QCD transition. This is 



in contrast to Ref. [14|. 
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FIG. 1. The sound speed (? s = (dp/dp) s during the QCD transition for the bag model (dotted 
line) and for the lattice QCD fit (full line). 
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FIG. 2. The modifications of the CDM density contrast ^4cdm = |#cdm|(I*/10) and of the 
radiation fluid amplitude ^4rad = (<5rad + ^rad) 1 ^ 2 due to the QCD transition (lattice QCD 
fit). Both quantities are normalized to the pure Harrison-Zel'dovich radiation amplitude. On the 
horizontal axis the wavenumber k is represented by the CDM mass contained in a sphere of radius 
ir/k. 
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FIG. 3. The evolution in conformal time r\ of the density contrast (5rad) and the peculiar 
velocity (~ V>rad) of the radiation fluid for the highest peak of Fig. 2 in uniform expansion (Hubble) 
gauge. During the QCD transition in the lattice QCD fit — marked by the 2 vertical lines — the 
velocity stays approximately constant and the density contrast grows linearly. The amplitude in 
the WKB regime is normalized to 1 long before the transition. 
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